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We demonstrate that it is possible to determine the coefficients of an all-order beta function linear 
in the anomalous dimensions using as data the two-loop coefficients together with the first one of the 
anomalous dimensions which are universal. The beta function allows to determine the anomalous 
dimension of the fermion masses at the infrared fixed point, and the resulting values compare well 
with the lattice determinations. 
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Understanding the non-perturbative dynamics of 
gauge theories of fundamental interactions constitutes 
a formidable challenge. Recently a considerable effort 
has been made to unveil the large distance conformal 
dynamics of these theories, see [1] for a recent review. 

The goal here is to prove the existence of an all-orders 
beta function similar, in shape, to the one provided by 
Ryttov and Sannino (RS) in [2|. 

We consider a generic vector-like gauge theory with N r 
Dirac fermions transforming according to distinct repre- 
sentations r of the underlying gauge group. The beta 
function of any gauge theory is a gauge independent 
quantity and therefore can only depend on gauge in- 
variant quantities. Besides, gauge invariance must be 
respected at each order in perturbation theory. Since the 
anomalous dimensions of the fermion masses are gauge 
independent the beta function can depend on them. We 
can therefore always write: 



P(a) = f(a,N r y r ,Yg), with r = l,...,p. 



where 



/3(a) 



da 
din [i 



and 



Yr 



dh\m r 
dlnp 



(1) 



(2) 



m r is the Dirac mass of each fermion species. y g is the 
gluon wave function renormalization anomalous dimen- 
sion satisfying the important relation: 



7s 



a 



(3) 



in the background field method ||3]-0. Note that the 
background field method does not fix uniquely the 
gauge. More precisely the quantities y r s and y g do 
not depend on the gauge fixing parameter while, on 
the the other hand the anomalous dimensions of the 
wave functions of the fermion fields do depend on this 
parameter. It will become soon clear that these anoma- 
lous dimensions are the minimal set of gauge invariant 
quantities needed to determine the beta function. 



We assume f(a,N r y r ,y g ) to be such that ji(a)/a 2 is 
linear in the anomalous dimensions: 



a 



a 
~2~j\ 



a + 



r 

Y^a r N r y r ■ 



(4) 



with r = l,...,p labeling matter transforming accord- 
ing to distinct representations of the underlying gauge 
group. We stress that the coefficients a, a r and a g are 
independent on a by assumption. We will prove below 
that, in perturbation theory, a scheme exists for which 
Eq. Q is valid. 
Using Q we find: 



a 



a a + Yl r=x a r N r yr 



2n 



(5) 



We will now show that it is possible to determine the 
p + 2 unknown coefficients a, a r and a g using the uni- 
versal coefficients of the two-loop beta function together 
with the universal coefficient of the anomalous dimen- 
sion of the mass for each representation. We henceforth 
introduce the two-loop beta function: 



2L 



a 



(6) 



with jSo an d f}\ the two universal coefficients: 



ft 



11 4fi 

yC 2 [G]--£T[r]IV,. / (7) 



-C 2 [G] 2 - fc 2 [G] £ T[r]N r - 2 £ C 2 [r]T[r]N r . 



(8) 



C2M is the quadratic Casimir and T[r] are the normal- 
izations of the generators in the representation r, and G 
indicates the adjoint representation. The explicit expres- 
sions for C2I/] and T[r] can be found in |1|. Expanding 
the all-orders beta function to two-loop and from the 
matching of the coefficients we find: 
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r 

a = jS , UgpO + ^ «r k r N r = jSj 



(9) 
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with k r = 3C 2 (r) a universal constant associated to the 
anomalous dimensions as follows: 



Yr 



(10) 



The remaining p + 1 coefficients, a g and a r in |4) can 
be determined from (9). We start by observing that, in 
general, the coefficients a g (N r ) and a r (N r ) can depend 
on the the number of fermion species N r . We start by 
considering the pure Yang-Mills case for which the only 
unknown coefficient a g is immediately deduced from (9| 
by setting N r = to be: 



ff M 17 



(11) 



and therefore the pure Yang-Mills beta function becomes: 



jgYM(a) = C 2 [G] 



(12) 



This is the RS beta for the pure Yang-Mills. The result- 
ing mass gap and the study of the poles of this theory 
have been analyzed recently in [6J. We also understand 
why the running of this beta function captures so well 
the lattice simulations performed using the Schrodinger 
functional scheme as shown in 1 2 1 . The reason is that the 
present beta function and the one derived on the lattice 
are both based on a background field method implying 
automatically the relation (3) in the two cases. 

To show how to derive the other coefficients we con- 
sider first the case of a single matter representation R for 
which the constraint in (9) reads: 



a g p +a R k R N R = jSi- 



(13) 



We therefore determine the coefficient a R , at the value 
N R for which asymptotic freedom is lost, i.e. j8o = 0, to 
be: 



a R (N R ) = 



k R N R 



(14) 



The two conditions 111 I and JT4} did not have, in prin 



ciple, to lead to a universal solution valid for any num- 
ber of matter fields. However, it is straightforward to 
show that a solution to |9) is obtained if a g and a R are 
constant, i.e. do not depend on N R , and equal to the 
values given by (TT) and (14) . In fact one can consider 
the two cases in which either a R does not depend on N R 
or the case in which a g is constant. Surprisingly these 
two limits lead to the same solution. The reason lies in 
the linearity of the two-loop beta function coefficients 
with the number of flavors N R . Thus the solution with 
a R and a„ constants independent on Nr is the most nat- 
ural choice. We note, however, that it is not the only 
logical possibility. In fact other solutions can be easily 
built, for example, by choosing any function a R (N R ) such 



that a R (N R ) = a R (0) = Pi(N R )/(k R N R ), and then choosing 
a g (N R ) to satisfy (9). The constraints on a R guarantees 

the solution to be regular at Nr = and Nr as we will 
discuss in more detail below. In the following we will 
only consider the solution with a R and a g constant and 
study its consequences. 

Given the two coefficients a g and a R the all-orders beta 
function j5}, for a single representation R, is determined 
to be 



a 



2tt 



1 



gYM 

a Pi 

2n p-™ 



(15) 



The expression above is, in shape, identical to the RS 
beta function [2J and, for Nr = 0, matches the Yang-Mills 

result. However for Nr — > Nr the RS does not reproduce 
the value of the anomalous dimension at the perturbative 
infrared fixed point. This occurs because the RS solution 
does not satisfy the condition (H) . We will show below 
that (14) is, in fact, a necessary condition if the perturba- 
tive expansion of the anomalous dimension at the fixed 
point has to be recovered. 



We note that, in general, if either 111 I or (14) are not 
satisfied, the coefficients a g or a R will have pole singu- 
larities at Nr = or Nr = Nr. In fact, for example, it 
is easily seen from |9) that if (H) is not satisfied then 
a g (N R ) ~ 1/jSo for Nr — » N R . This generalizes to any 
number of matter representations. 

Note also that for the solution proposed here the de- 
nominator does not depend on the number of flavors. 
Finally the value of the anomalous dimension at the 
fixed point is corrected with respect to the RS result. 
We rewrite the beta function after having evaluated the 
different coefficients: 



m 

a 



11C 2 [G]-2T[R]Nr(2 + Ar7r) 



6tt 



1 " 27?TT C 2[ G ] 



with 



1 + 



7 C 2 [G] 
11 C 2 [R] 



(16) 



(17) 



Interestingly in [7J the same form of the beta function 
was also suggested among a one-parameter family of 
solutions of (9), partially motivated from holography 
and assuming the RS shape of the beta function. It is 
worth stressing that our main point here is to understand 
the assumptions behind a field theoretical justification 
for the RS form of the beta function and the uniqueness 
of the solution introduced here. 

The new analytical expression of the anomalous di- 
mension of the mass at the infrared stable fixed point is 
obtained by setting the beta-function to zero and reads: 



11C 2 [G] - 4T[R]N f 



Po(N R ) N Rl 
Yr = ~ — — — k R 



(18) 
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The RS result is instead: 



, RS _ 11C 2 [G]-4T[R]N R 
7r 2N r T[R] 



(19) 



It is useful to compare the two anomalous dimensions: 



Yr 



11 



11 + 7 



C 2 [G] 
C 2 [R] 



(20) 



This shows that the corrected anomalous dimension at 
the fixed point is smaller than the one predicted ear- 
lier, which agreed roughly with the Schwinger-Dyson 
results. 

We now show that if a r differs from fl4) the anoma- 
lous dimension at the fixed point obtained via ||4} does 
not reproduce the perturbative value. The latter is ob- 
tained using the two-loop beta function and the one-loop 
gamma function which gives: 



Yr 



perturbative = + Q ^ 2) 



(21) 



On the other hand from the proposed beta function we 
obtain (Tl8>: 



^-(a r (N R )N R )~ 1 p +O(li 2 ). (22) 



YR 



Requiring, for consistency, equations pT) and | [22) to 
agree in perturbation theory allows to recover exactly 
(111. This new condition implies that only the regular 
solutions, i.e. non singular in jSo, are acceptable when 
solving 

Next, we would like to show that it is always possible 
to find a renormalization scheme in which | |T5| | holds. To 
this goal we introduce the general perturbative expan- 
sion for the anomalous dimensions and the beta function: 



n=0 



Yr 



CO 

n=l 



(23) 
(24) 



and by inserting the equations above in |[T5j we derive 
the infinite set of relations: 



fin ~ f>n-\ dg 

Yn-i = — r; , n > 1 , 

N R a R 



(25) 



with yo = k R - 3C 2 (R). The transformation laws be- 
tween the beta function and the anomalous dimensions 
above and the ones in another mass-independent renor- 



malization scheme /3(a) and f R (a) are: 



h„a n , with hi = 1 , 



(26) 



«=i 



in 



= m z m (a) = m ^ 4«" , with 4 = 1, (27) 



Yr(oc) = f R (a) - 



n=0 



<91nz„ 
da 



(28) 
(29) 



yo, |3n and jSi are unchanged being universal. The other 
coefficients do depend on the transformation and we 
report the explicit form of the first coefficients: 



h 
Y\ 



4n 2 p (h 2 2 - h 3 ) + 2nh 2 ^ + /3 2 , 
27rfijS + 2nh 2 y + fi ■ 



(30) 
(31) 



Imposing the relation (25J between y\, $\ and f> 2 yields: 

4ti 2 j So(^ - h) - 2na R N R { 1 p - a g ^ + 2nh 2 f, x + 

+p 2 - 2na R h 2 N R y - UrNr^ = . (32) 

To each successive order in perturbation theory, two 
more coefficients (one from a and one from z m ) appear 
and only one relation coming from ( |2"5"|l s hould be im- 
posed. Therefore we can always write \V5\ as a series in 
perturbation theory. 

The new anomalous dimensions at the fixed point in 
|T8"l l are in better agreement with the lattice determina- 
tions, which tend to be smaller than the RS estimate. The 
value of the anomalous dimensions at the fixed point is 
a physical quantity and therefore does not depend on 
the renormalization scheme. This fact can be directly 
verified using Eq. ( (29) . We quote below the value ob- 
tained for different theories on the lattice which have to 
be taken cum granum salis given that they are still subject 
to large systematic errors difficult to estimate. 

The anomalous dimension of the mass for the Min- 
imal Walking Technicolor theory |8j O corresponding 
to an SU(2) gauge theory with two Dirac flavors in the 
adjoint representation is now predicted to be 7mwt = 
11/24 0.458 rather than 0.75 as predicted earlier. 
This result compares well with the latest lattice result 
Ymwi ~ 0-49(13) [10 1, while a more conservative lattice 
estimate indicates the fixed value to be smaller than 0.56 
1111 . The anomalous dimension of the mass for the Next 
to Minimal Walking Technicolor theory (NMWT) cor- 
responding to an SLT(3) gauge theory with fermions in 
the two index symmetric representation is predicted to 
be 7nmwt = 143/173 - 0.827 which is closer to the one 
obtained via first principle lattice simulations in 1 12J. 

For 3 colors and 10 flavors in the fundamental repre- 
sentation we find yFund - 0.53 versus 1.3 obtained via 
the RS result. 

To provide further support to our results we compare 
the prediction for the anomalous dimension at the fixed 
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tion: 



FIG. 1: Comparison of the predicted value of anomalous di- 
mension at the fixed point from the proposed all-order jS func- 
tion and the one coming from perturbation theory up to 4 
loops. As a representative gauge theory we used here the SU(2) 
gauge group with fermions in the adjoint representation. Sim- 
ilar curves are obtained for other gauge theories. From the top 
curve down: 2, 3 and 4-loop perturbative results. Dashed line: 
anomalous dimension coming from the all-order /5 function. 



point given in Eq. with the results of perturbation 
theory, known up to 4 loops. We find a remarkable 
agreement between the two results. This is a surprising 
agreement given that our proposal makes use only of the 
2-loop universal coefficients of the |3 and y functions. To 
better appreciate the striking agreement with perturba- 
tion theory, we plot in Fig. [lithe fixed point value of )Ar 
from Eq. ( jl"8) and from perturbation theory at 2, 3 and 
4-loop, in the case of SU(2) gauge theory with adjoint 
fermions. The prediction coming from the proposed f>- 
f unction nicely envelops the perturbative results. Simi- 
lar results are obtained for other gauge theories. We tab- 
ulate in the Appendix the values of y R corresponding to 
the gauge groups SU(2), SU(3) and SLT(4) with fermions 
in the fundamental and the two indices representations. 

Using the linearity in the number of flavors of the 
first two-coefficients of the beta function it is straightfor- 
ward to generalize the beta function to any number p of 
fermionic species, the result being: 



a 



a llC 2 [G]-2Z P r=1 Tlr]N r (2 + A r y r 



6tt 



a_17 
In 11 



C a [G] 



(33) 



with A,- given by ( [17} . Our results generalize immedi- 
ately to the supersymmetric case [ 1 3 1 by simply replacing 
the corresponding quantities with the supersymmetric 
ones in ( 11 1 and 1 14 1 and obtain, for a single representa- 



j3 s (a) = a gg + 2T[R]N R y R 
a ~ 2n 1-§C 2 [G] 



(34) 



where the anomalous dimension of the mass is defined 
with an overall sign difference with respect to the non- 
supersymmetric case and the S apex means supersym- 
metric, for example, /3g = 3C 2 (G) - 2T[R]N R . We find re- 
markable that the beta function is shape invariant when 
going from the nonsupersymmetric to the supersymmet- 
ric parent theory. This was not the case for the RS beta 
function. 

We have argued that the beta function for any non- 
supersymmetric vector like gauge theory with fermionic 
matter as well as its supersymmetric version can be 
written in the universal form (15) . The latter is ob- 
tained solely in terms of the universal coefficients of the 
two-loop beta function and the universal coefficient of 
the first term in the coupling constant expansion of the 
anomalous dimension of the mass operator. 



From 1 15 1 it is possible to read the values of the anoma- 
lous dimension of the mass inside the conformal win- 
dow, which is found to be in fair agreement with the 
numerical estimates obtained via numerical simulations. 

We thank Oleg Antipin, Marco Nardecchia, Thomas 
A. Ryttov and Joseph Schechter for helpful comments. 



Appendix A: Perturbative versus all orders results 

In this section we report the values of the fixed point 
anomalous dimension of the fermion mass for different 
gauge theories as obtained from Eq. (18) y* AO , and from 
perturbation theory at 2-loop yl, 3-loop y* 3 and 4-loop 
y* v We also report the scheme dependent value of the 
coupling constant A* = a* /An corresponding to the zero 
of the perturbative /3 function at 2, 3 and 4-loop. The 

perturbative expressions are obtained in the MS scheme 
IIT4T IT6 I. The gauge theories considered in this Appendix 
are based on gauge group SU(2), SU(3) or Slf(4) with 
fermions in the fundamental, adjoint, 2-index symmetric 
or anti-symmetric representation. 

It is clear that for sufficiently low number of flavors, 
as expected, the four-loop approximation breaks down. 
This is the case, for example of the 7 and 6 flavors for 
the SU(2) with fundamental fermions where one ob- 
serves wild fluctuations of the fixed point couplings and 
anomalous dimensions at different orders. Notably the 
all-orders results are, on the other hand, stable. 
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n f 
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TABLE II: Same as Table [I] for the SU(3) gauge theory with 
fundamental fermions. 
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SU(4) Fundamental 



n f 


iif/rif 


A\ 


y\ 






y\ 


K 


y\ 




11.000 





.500 


1 


.110 241.000 


.077 


1 


.830 





.075 


-0.227 0. 


.848 


11.500 


0. 


.523 





.463 


44.000 





.068 


1 


,400 





.067 


-0.053 0. 


.775 


12.000 


0. 


.545 





.282 


17.300 





.060 


1 


.110 





.060 


0.058 


0. 


.707 


12.500 


0. 


.568 





.197 


8.980 





.054 





.904 





.055 


0.136 


0. 


.645 


13.000 


0. 


.591 





.147 


5.380 





.048 





.755 





.050 


0.192 


0. 


.587 


13.500 





.614 





.115 


3.520 





.043 





.642 





.046 


0.232 





.534 


14.000 





.636 





.092 


2.450 





.039 





.552 





.041 


0.259 


0. 


.485 


14.500 


0. 


.659 





.075 


1.770 





.035 





.480 





.038 


0.275 


0. 


.439 


15.000 


0. 


.682 





.062 


1.320 





.032 





.420 





.034 


0.281 


0. 


.396 


15.500 


0. 


.705 





.052 


1.010 





.028 





.369 





.031 


0.278 


0. 


.356 


16.000 


0. 


.727 





.043 


0.778 





.025 





.325 





.027 


0.269 





.318 


16.500 


0. 


.750 





.036 


0.610 





.023 





.286 





.024 


0.254 





.283 


17.000 


0. 


.773 





.031 


0.481 





.020 


0. 


.251 





.022 


0.234 


0. 


.250 


17.500 


0. 


.795 





.026 


0.380 





.018 





.218 





.019 


0.211 


0. 


.218 


18.000 


0. 


.818 





.021 


0.301 





.015 





.189 





.016 


0.187 


0. 


.189 


18.500 





.841 





.017 


0.236 





.013 





.161 





.014 


0.162 


0. 


.160 


19.000 





.864 





.014 


0.183 





.011 





.134 





.012 


0.136 





.134 


19.500 


0. 


.886 





.011 


0.139 





.009 





.109 





.009 


0.111 


0. 


.109 


20.000 


0. 


.909 





.008 


0.102 





.007 





.085 





.007 


0.086 


0. 


.085 


20.500 


0. 


.932 





.006 


0.071 





.005 





.063 





.005 


0.063 





.062 


21.000 


0. 


.955 





.004 


0.044 





.004 





.041 





.004 


0.041 


0. 


.040 


21.500 


0. 


.977 





.002 


0.021 





.002 





.020 





.002 


0.020 





.020 


22.000 


1. 


.000 





.000 


0.000 





.000 





.000 





.000 


0.000 


0. 


.000 



TABLE III: Same as Table [I] for the SU(4) gauge theory with 
fundamental fermions. 



SU(2) Adjoint 

n f n f /n f A\ y' 2 A" 3 y\ A* y\ y\ 
1.500 0.545 0.179 5.370 0.083 1.920 0.068 1.300 1.020 
2.000 0.727 0.050 0.820 0.037 0.543 0.036 0.500 0.458 
2.500 0.909 0.011 0.139 0.010 0.127 0.010 0.127 0.122 

SU(3) Adjoint 

1.500 0.545 0.119 5.370 0.055 1.920 0.049 1.520 1.020 
2.000 0.727 0.033 0.820 0.024 0.543 0.025 0.523 0.458 
2.500 0.909 0.007 0.139 0.007 0.127 0.007 0.128 0.122 



SU(4) Adjoint 



1.500 0.545 0.089 5.370 0.041 1.920 0.038 1.620 1.020 
2.000 0.727 0.025 0.820 0.018 0.543 0.019 0.532 0.458 
2.500 0.909 0.005 0.139 0.005 0.127 0.005 0.128 0.122 

TABLE IV: Same as Table|I]for the gauge theories SU(2), SU(3) 
and SU(4) with adjoint matter. 



SU(3) 2-index Symmetric 



n f n f ln f A\ y\ A' 3 y' 3 A\ y\ y\ Q 

1.500 0.455 0.261 23.800 0.081 5.290 0.061 3.590 1.530 

2.000 0.606 0.067 2.440 0.040 1.280 0.037 1.120 0.827 

2.500 0.758 0.025 0.639 0.020 0.474 0.020 0.466 0.407 

3.000 0.909 0.007 0.144 0.006 0.133 0.006 0.133 0.127 



SU(4) 2-index Symmetric 



1.500 0.409 0.441 111.000 0.079 10.400 0.056 6.750 1.850 

2.000 0.545 0.077 4.820 0.039 2.080 0.035 1.790 1.060 

2.500 0.682 0.030 1.210 0.021 0.776 0.021 0.746 0.596 

3.000 0.818 0.012 0.381 0.010 0.313 0.010 0.315 0.284 

3.500 0.955 0.002 0.064 0.002 0.062 0.002 0.062 0.061 

TABLE V: Same as Table [I] for the gauge theories SU(3) and 
SU(4) with 2-index symmetric matter. 
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SU(4) 2-index Antisymmetric 



n f 


n-f/rif 


A* 


y\ 






rl 




y* 


7'ao 


5.000 


0.455 


4 


.000 4200.000 


.089 


4 


.150 





,130 


-4.940 


1.190 


5.500 


0.500 





.361 


38.000 





.067 


2 


.210 





,079 


-0.026 0.991 


6.000 


0.545 





.172 


9.780 





.053 


1 


.380 





,061 


0.293 


0.826 


6.500 


0.591 





.105 


4.170 





.043 





.952 





,049 


0.402 


0.686 


7.000 


0.636 





.071 


2.190 





.035 





.695 





,040 


0.435 


0.566 


7.500 


0.682 





.050 


1.290 





.028 





.525 





,032 


0.418 


0.462 


8.000 


0.727 





.036 


0.802 





.023 





.402 





,025 


0.368 


0.372 


8.500 


0.773 





.026 


0.515 





.018 





.306 





,020 


0.302 


0.291 


9.000 


0.818 





.018 


0.331 





.014 





.228 





,015 


0.232 


0.220 


9.500 


0.864 





.012 


0.206 





.010 





.161 





,010 


0.164 


0.156 


10.000 


0.909 





.007 


0.117 





.006 





.101 





,007 


0.103 


0.099 


10.500 


0.955 





.003 


0.051 





.003 





.048 





,003 


0.048 


0.047 


11.000 


1.000 





.000 


0.000 





.000 





.000 


0. 


,000 


0.000 


0.000 



TABLE VI: Same as Table [rj for the SU(4) gauge theory with 
2-index antisymmetric matter. 



